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We investigate the spectral shape of a high-energy diffuse photon emitted by evaporating primor- 
dial black holes (PBHs) in the Randall-Sundrum type II (RS2) braneworld. In their braneworld 
scenario, the nature of small PBHs is drastically modified from the ordinary four-dimensional case 
for the following two reasons, (i) dropping Hawking temperature, which equivalently lengthens the 
lifetime of the individual PBH due to the change of space-time topology and (ii) the effective increase 
of the total amount of PBHs caused by accretion during the earliest part of the radiation-dominated 
epoch, the brane high-energy phase. From studies of the expected spectral shape and its dependence 
on braneworld parameters, we obtain two qualitatively distinctive possibilities of constraints on the 
braneworld PBHs from the observations of diffuse high-energy photon background. If the efficiency 
of accretion in the high-energy phase exceeds a critical value, the existence of the extra dimension 
gives a more stringent upper bound on the abundance of PBHs than the 4D case and a small length 
scale for the extra dimension is favored. On the contrary, in the case below the critical accretion 
efficiency, we find that the constraint on the PBH abundance can be relaxed by a few orders of 
magnitude in exchange for the existence of the large extra dimension; its size may be even bounded 
in the region above 10^^ times 4D Planck length scale provided the rest mass energy density of the 
PBHs relative to energy density of radiation is actually larger than 10~^^ (4D upper bound) at their 
formation time. The above analytical studies are also confirmed numerically, and an allowed region 
for braneworld parameters and PBH abundance is clearly obtained. 

PACS numbers; 98.80.Cq, 04.50.-|-h, 04.70.Dy 



I. INTRODUCTION 

It is widely believed that the local density contrast 
in the inflationary early Universe caused gravitational 
collapse, which resulted in the subsequent formation of 
plenty of mini black holes, the so-called primordial black 
holes (PBHs) Him Such small black holes are 

thought to be only able to have formed in the early Uni- 
verse. There are even other mechanisms that could cause 
the PBH formation 6, 7, 8]. 

Since Hawking discovered the possibility that black 
holes can emit particles while decreasing their masses 
0, 01 , investigating the energy spectra of emission from 
PBHs has been one of the vital ways to obtain informa- 
tions on the earliest era of the Universe. Small black 
holes in the Universe radiate like blackbodies at every 
moment in their lives and leave traces in their relics. In 
particular PBHs which have lifetimes around the present 
cosmological age have been emitting only massless par- 
ticles, whose optical depths are much less than unity; 
they are in principle directly detectable. Hence compar- 
ing predicted spectra with present observations, we can 
obtain information on the early Universe via the abun- 
dance of PBHs. One of the earliest works that applied 
the theory of black hole evaporation into astrophysics 



was carried out by Page and Hawking |ll|, who consid- 
ered the diffuse 7-ray background to probe the number 
density of PBHs that have evaporated just now. With the 
same motivation, a number of ar gum ents on the amount 
of PBHs have been reported [llllllT3|. 

More than two decades after Hawking's discovery, it 
was suggested by Randall and Sundrum 15] that our 
Universe is a 4D Lorentz metric hypersurface, called 
the brane, embedded in Z2 symmetric 5D anti-de Sit- 
ter (AdS) spacetime, called the bulk. The so-called RS2 
braneworld model, which we concentrate on in this paper, 
is one of the probable candidate models for our Universe 
they suggested. It assumes that the number of branes 
is only one and the brane has positive tension. In their 
scenario, ordinary matter is confined on the brane while 
4D gravity is reproduced due to effective compactifica- 
tion. The bulk has a typical length scale that originates 
from AdS curvature radius, /. This length scale draws a 
short distance boundary on the brane; deviation from or- 
dinary gravity appears below it, while 4D Newton law is 
recovered as far as we consider a sufficiently larger scale 
than I 0, ^1 . The present table-top upper bound on 
the bulk AdS radius I was obtained by measuring short- 
distance gravitational force as / < 0.2 mm 17] . In other 
words, there is a possibility that brane signature is di- 
rectly seen by observing submillimeter scale gravitational 
interaction. 
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From such a point of view, a black hole, the most 
strongly gravitationally bounded object, is one of the 
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most appropriate probes for submillimeter physics. Ac- 
tually, in the framework of the braneworld it has been 
revealed that the nature of black holes with small masses 
highly differs from the ordinary 4D one. In particular 
small black holes with Schwarzschild radius rs ^ I be- 
have as essentially 5D objects and their mass-radius re- 
lation is modified [11 El IM El 113 ■ On the other hand 
there have been discussions on the black hole formation 
in respect of the particle scattering processes [2I |23| . 
but we naturally expect that even such small black holes 
were produced by gravitational collapse in the early uni- 
verse I2I. Thus we have to revisit the particle emission 
process of primordial black holes taking into account how 
non-4D character modifies it. In addition, the RS2 sce- 
nario implies a nonstandard high-energy reginie before 
the standard radiation-dominated era |3 HI HI 113 
(see, e.g., 28] and references therein). It was recently 
argued by Guedens, Clancy, and Liddle and Majumdar 
[29. 30, 3^ that through the brane high-energy phase 
accretion of surrounding radiation onto PBHs is so effi- 
cient that they grow in their masses. On the other hand, 
such accretion can be neglected as the Universe enters 
the standard radiation-dominated era. 

Furthermore, taking into account the above five- 
dimensional signatures, Clancy, Guedens, and Liddle [s^l 
discovered that the shape of diffuse photon spectrum 
emitted from PBHs is drastically modified from that cal- 
culated by Page and Hawking in the 4D case, especially 
the peak energy can drop by three orders of magnitude. 
Therefore it was revealed that we can obtain crucial in- 
formation on the earliest universe in the braneworld sce- 
narios by investigating those spectral shapes carefully. 
However, what was obtained by Clancy et al. is only 
a rough shape of the photon spectrum. Although they 
even estimated an upper bound on initial mass fraction 
of PBHs, it was only for one specific parameter set where 
the effect of accretion was ignored. In other words they 
argued about a revised upper bound on the PBH abun- 
dance standing on the theory of 5D evaporating black 
holes, but 5D brane cosmology was not significantly taken 
into consideration. 

On such a background our purpose in this paper is to 
reveal the detailed dependence of spectral shape of emit- 
ted photons on the size of extra dimension I and accre- 
tion efficiency F. For this purpose we first derive how the 
PBH mass spectrum is modified by accretion, which has 
not been carefully considered before. Next we estimate 
the emitted photon spectrum taking the modified PBH 
mass spectrum into account, and compare them with 
the diffuse high-energy photon background spectrum ob- 
served by HEAO-1 A2 33], HEAO-1 A4 Q, SMM [H, 
COMPTEL [13, and EGRET In the analysis, we 

particularly pay attention to the fact that the observed 
spectrum obeys nearly a power-law of index — 1 in almost 
all energy regions. It turns out later that this index de- 
termines a critical accretion efficiency. Additionally we 
also concentrate on the origin of difference between the 
new 5D constraints and the old 4D one. We find that 



there is a kind of critical size around Z « 5 x lO^^l^ in the 
framework of evaporating PBHs. 

For further analysis we even numerically study dif- 
fuse high-energy photon spectra for various possible 
braneworld parameter sets. We use throughout this 
paper the newest observational cosmological parame- 
ters based on Wilkinson Microwave Anisotropy Probe 
(WMAP) observations 37], which confirmed the pre- 
ceding suggestions of a A-dominated universe obtained 
from measurements on high-redshift type la supernovae 
m [33. We also take into consideration precisely the 
effect of a mass decrease of PBHs in the evaporating pro- 
cess. Hence expected shapes of diffuse photon spectrum 
from braneworld PBHs are obviously shown. 

The contents of this paper are as follows. First we 
briefly review in Sec. HTlcosmologv and behavior of evap- 
orating black holes in the RS2 braneworld scenario. In 
Sec, mil formulations to describe PBHs in the brane uni- 
verse are shown. Those are basically according to the 
preceding papers [H [sl |^ [sj. After the theoreti- 
cal preparations, we analyze general dependences of the 
PBH mass spectrum and resulting photon spectrum on 
the braneworld parameters in Sec. IIVI In this section, 
our main purpose, new constraints on braneworld pa- 
rameters and their implications are obtained. Finally in 
Sec. we show numerical results on the calculation of 
diffuse photon spectra for a wide range of braneworld 
parameters and refine and visualize our analysis in the 
previous section. Conclusions and discussions are pre- 
sented in Sec. IVII 



II. COSMOLOGY AND BLACK HOLES IN THE 
BRANEWORLD 

A. Brane cosmology 

Hereafter we basically use 4D natural units, that is 
h = c = kB = 1 throughout this paper. On the other 
hand 4D Planck scale quantities defined as M4, 14, ti, and 
T4 usually appear in expressions because the braneworld 
scenario, which we consider here, is a higher-dimensional 
unified theory. 

We introduce a five-dimensional fundamental scale K5 
of mass dimension —3/2, which appears in the 5D Ein- 
stein equation as 

r<(5) _ A „(5) , ^2rp(5) /r,-,N 

'^MN— ^^^y M N ^ '^b-'- M N ■ \'^-^) 

Here 5^^^^ is 5D metric, G^^jj^ and r^]^ are the 5D Ein- 
stein tensor and energy-momentum tensor, respectively. 
A5 is the negative cosmological constant in the bulk. 
Expansion-law on the brane is determined by this five- 
dimensional equation with an appropriate boundary con- 
dition on the brane, so-called Israel's junction condition. 
In this framework, brane cosmological expansion is mod- 
ified from the 4D case due to the existence of typical bulk 
length scale, the AdS curvature radius /, and energy scale 
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on the brane, the brane tension A lia|2J,|25|,|2a|27|- The 

AdS curvature radius is defined by bulk energy scales as 



l = J- 



6 



(2.2) 



The quantity I expresses how much the bulk has physi- 
cal volume, hence according to the size of I we mention 
whether the extra dimension is large or small. 

Using the 5D Einstein equation in the bulk and the 
boundary condition on the brane under Z2 symmetry, 
one can obtain the 4D Friedmann equation on the brane 
in a modified form: 



Stt 
3M| 



(-1;) 



Pkk 



A4 
3 



K 
7^' 



(2.3) 



where a is the scale factor, H is the Hubble constant, p 
and Pkk are the energy density of ordinary matter field 
and of the dark radiation, respectively, and K is the spa- 
tial curvature of the brane. M4 and A4 are the induced 
Planck mass and cosmological constant on the brane re- 
lated to 5D quantities as 



A4- ^ 
2 



A5 + 
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It is needed for these quantities to match observed 4D 
values. The energy conservation law is the same form as 
the 4D one: 



p + 3H{p + p) = 



(2.5) 



(see, e.g., [2a| and references therein for a review). 

Here we evaluate the terms in the above modified 
Friedmann equation. As we are interested in a case which 
is connected with the well-established standard 4D uni- 
verse, each term is bounded or even determined by cur- 
rent observations. As for the dark radiation term, it was 
argued that this quantity must be sufficiently smaller 
than ordinary radiation at the time of nucleosynthesis 
[27! 5^ , thus we drop this term. Next the observations 
of WMAP imply that spatial curvature of the brane is 
small [33, hence we also set ivT = 0. Though WMAP 
also revealed that the 4D cosmological constant A4 takes 
a significant role in the present Universe, it can be ig- 
nored in the early Universe. We require A4 to vanish 
except when the present epoch of the Universe is con- 
sidered. Therefore one finds formal Friedmann equation 
(|2.3|) results in the following form: 



H^^^p(l + J^)+^ 
ZMVK 2\) 3 



(2.6) 



The above expression differs from the ordinary 4D one 
with respect to the existence of a p-square term. This 
term is only relevant when energy density is sufficiently 
higher than the brane tension, namely in earliest era of 
the Universe. Then we can split the cosmological evo- 
lution in the early Universe into two phases; the ear- 
lier one is the nonstandard high energy regime in which 



evolution of the scale factor a{t) obeys a{t) oc t^/^, and 
the latter is the standard radiation-dominated universe, 
where a{t) oc t^/^. The transition time between them is 
tc = 1/2. 

Physical quantities governing cosmological expansion 
should be determined by observation. Thus the proce- 
dure of determining the profile of scale factor a{t) is as fol- 
lows. To begin with, we normalize the scale factor so that 
its present value becomes unity: a(to) = 1- Then we inte- 
grate the Friedmann equation to the past using WMAP 
data, such as to = 13.7 Gyr, h = 0.71, ilm^o^^ — 0.135, 
^^A,o = 0.73, and Zeq = 3233. On the other hand, as for 
radiation energy density, observed data taken by Cosmic 
Background Explorer (COBE) 0| is used as the present 



value: ilifih 



4.15 X 10" 



Then all we have to do 



is go back to the time when ^ holds and con- 
nects to the Einstein-de Sitter universe. The choice of 
"connection time," t^, is artificial, and hence eventually 
may cause an error in other quantities; we carefully chose 
an appropriate t\ so that the time of matter-radiation 
equality ^oq, which is determined by a^q with ^a, has no 
meaningful numerical error. As a consequence the profile 
of scale factor is determined as 



a{t) 



1/2 
/2/3 



for i < ic 



(2.7a) 



for tc < i < icq (2.7b) 



for teq<t< tA, (2.7c) 



where the value of icq is calculated as 72.6 kyr and a 
prime denotes redefinition of time origin {f = t[,^ when 
t = ^eq)- Note that although we often make use of 
these asymptotic forms there is even an analytic solution 
through the whole radiation-dominated era: 



i{t) 



1/2 
oq 



(2.8) 



It is also utilized to describe near the transition time t^. 

Using the above profile of a{t), other cosmological 
quantities, energy density p{t) , Hubble radius i?H [t) , and 
horizon mass M-H_{t) are calculated. We are particularly 
interested in their expressions in the early Universe: 



Pit) 


3Mi 
327rtci' 


Ruit) = At, 


Mait) 


= 8M^- 


for t < tc, 


Pit) 




Ruit) = 2t, 




= Mlt 


for t > tc- 



(2.9a) 



(2.9b) 
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B. Brane black holes 

In the braneworld scenario, a massive object localized 
on the brane with typical radius r <^ I should see the 
background space-time as effectively flat 5D Minkowski 
space [i^ |2y, |2l|, |23 • Therefore a sufficiently small black 
hole is no longer an ordinary 4D one, but can be approxi- 
mately described as 5D Schwarzschild solution. In terms 
of genera l sp herical coordinate, the 5D Schwarzschild 
metric is 143 



ds^ 



with 



-C{r)At^ 



C{r) 



1 



-dr^ 



C{r\ 



r 



(2.10) 



(2.11) 



where rs is the Schwarzschild radius and dfi d is the line 
element on unit ZJ-sphere. The Schwarzschild radius in 
the braneworld is given as 



III. FORMULATIONS 
A. PBH formation 

It has been argued that in the early Universe PBHs 
are formed due to gravitational collapse caused by den- 
sity perturbation. Though there are theoretical difficul- 
ties on the process of gravitational collapse,^ we here 
simplify the situation; when a superhorizon scale density 
perturbation with adequate amplitude enters the Hub- 
ble horizon at t = ti, it instantly forms a black hole with 
some fraction of horizon mass. Hence the formation mass 
M{ is given by 



Mi = /MH(ti), 



(3.1) 



where / expresses the fraction of horizon mass collaps- 
ing into the black hole. It can be confirmed that / is 
order unity by virtue of Jeans analysis and that such a 
black hole has a Schwarzschild radius nearly covering the 
Hubble horizon 12911. 




(2.12) 



The induced 4D metric near the event horizon is obtained 
by projecting the bulk 5D metric onto brane, 



dsl 



-C(r)di2 + ^dr^ + r'^dnl. (2.13) 



On the other hand general D-dimensional Hawking tem- 
perature of a Schwarzschild black hole is given in terms 
of its radius as l20l 1431 



D 



47rrs 



(2.14) 



Combining these expressions, one finds the modified mass 
to Hawking temperature relation. 



BH 



3 / J_ 
327r W4 



-1/2 



M 
M4 



-1/2 



T4 



While ordinary 4D expressions are 
C(r)(4D) = 1-^, 



rs(4D) 
Tbh(4D) 



M4 

-f-V 



(2.15) 

(2.16) 
(2.17) 

(2.18) 



As far as we consider such small black holes for which 
the above 5D approximations are valid, the temperature 
of a black hole with some fixed mass M is to be in general 
lower than in the 4D universe. We see in the next section 
how this change affects the spectra of particles emitted 
from black holes. 



B. Accretion on PBHs 

Radiation field localized on the brane sees a black hole 
approximately as a disk-shaped object with radius rs^cff, 
where rs^eff is the effective D-dimensional Schwarzschild 
radius given by [44| 



D - 1 



l/(D-3) 



D - 1 
L» - 3 



1/2 



rs- 



(3.2) 



Consider an infinitesimal time interval [t, t+dt] . Then the 
corresponding gain of black hole mass due to absorption 
of surrounding radiation is described as [s^, |^ 



dM 
"dT 



(3.3) 



where p{t) is energy density of radiation. Here we in- 
troduced a parameter _F, which determines accretion ef- 
ficiency. If the geometrical optics approximation holds 
well, namely an incidental relativistic particle can be re- 
garded as a collisionless point particle and its spin can be 
ignored, F should be almost unity, but otherwise smaller. 

Substituting relevant expressions into the above differ- 
ential equation and integrating it, we obtain the mass 
evolution relation in the brane high energy phase as 

MM 



M{t) 



2F/TT 



Mi for t<tc. 



(3.4) 



^ Actually there is a paper showing critical phenomena in the for- 
mation of PBHs 43]. However, it is still reasonable to take 
/ = 0(1) in Eq. i5H . 
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We define by the above formula Mc = M{tc). 

The effect of evaporation was ignored in the above 
derivation. Thus one may suspect the above expression 
is not vahd when accretion does not work, say F — 0. 
However, evaporation in such a short period naturally 
has little significance for long-life PBHs which we will 
consider later; we can safely make use of the above ex- 
pression regardless of the value of F. 



C. Evaporation 

The differential particle emission rate from a D- 
dimensional static noncharged black hole with mass M 
is given by Hawking's formula: 



AN, 
dt ^' 



a,{E) 



dk 



D-l 



eME/TBn) ± 1 (27r 



i-D-l ■ 



(3.5) 



where the subscript j indicates particle species, gj is the 
number of internal degrees of freedom of the particle, 

E ~ + mj is the relativistic energy including rest 

mass rrij, aj is the absorption cross section, the sign in 
the denominator expresses which statistics the particle 
obeys, and Tbh is the Hawking temperature related to 
the black hole Schwarzschild radius given in Eq. H2.14|l . 
Using the above formula under relativistic limit, the total 
mass decreasing rate of a spherically symmetric black 
hole results in 



dM 
'df 



E 
E 



E 



E- — Me 



dEdt 



9j 



9j 



nD-2E'^'^dE 



eME/TBn) ± 1 (2^)^-1 
exp{E/TBii) ± 1 (2^)^-1 



cS.dTbyi 



(3.6) 



Here we introduced some quantities following l29|; go is 
an effective degree of freedom given by 

2°-^ - 1 

go .9-D, boson H — 1 .^D.fcrmionj (3.7) 

ao is the D-dimensional Stefan-Boltzmann constant 



4(27r)^ 



-r(^)C(^), 



(3.8) 



and AcH^D is the {D — 2)-dimensional effective surface 
area 

Ac«M = ^D-2rg~^ jj, (3.9) 

where flD-2 is the area of the (D — 2)-dimensional unit 
sphere, and r(Z?) and C(^) gamma and zeta func- 
tions, respectively. 



To investigate the evaporation process of black holes 
in the braneworld, one should consider degrees of free- 
dom on the brane and in the bulk separately because 
the former takes 3D phase space and 2D black hole sur- 
face area while the latter takes 4D phase space and 3D 
surface area. Thus the mass shedding formula takes the 
form below: 

= — ffbranc^cff,4rBjj — 5bulk^cff,5'7BH • (3.10) 

Note that we must take Acs^i — Trrg ^ since 
the Schwarzschild radius is determined by the five- 
dimensional metric. 

From the above arguments, we obtained the fundamen- 
tal differential equation which describes mass-loss rate of 
black holes. Thus by integrating it, we know how the 
mass of a black hole decreases and therefore how the 
Hawking temperature rises during its lifetime. (A pos- 
sible deviation from the above formula is discussed in 

One theoretical uncertainty is on the treatment of ef- 
fective degrees of freedom for high temperature black 
holes with Tbh ^ "Itt, where m,^ ~ 135 MeV is the 
mass of pion as the lightest hadron. However, the heat- 
up rate of a black hole in its early phase is so slow and 
once Tbh gets as high as ttIt^, a black hole "blows up" 
and ends its life almost instantly. Thus in the following 
analysis, we can assume an evaporating black hole holds 
approximately the same temperature as the initial value 
through its lifetime, and hence the effective total degrees 
of freedom, gcs, is constant. With low temperature as- 
sumption (Tbh ^ 100 MeV), we can conclude that the 
relevant particle species to calculate gcs for those PBHs 
surviving until present epoch are massless species and at 
most lightest (anti-)leptons, i.e., photon, neutrinos, and 
electron/positron. In the last part of this section, we will 
check the consistency of the low temperature assumption 
and, for concreteness, estimate the fraction of photon en- 
ergy emitted after the assumption breaks down. On the 
other hand we still do not have enough knowledge to pre- 
cisely determine the spin-dependent characters of parti- 
cle emission. We here take the values obtained in the 4D 
case by Page ^ consequence, the mass shedding 

rate is given as 



dt 2 yhj V Mi J ti ' 



(3.11) 



We here cite the value of the effective degree of freedom 
goff in 29]: 

ffeff ~ 0.023 for only massless 

~ 0.032 for massless and electron/positron. 

Difference between them is almost negligible. The solu- 
tion for Eq. H3.11|) is easily obtained as 



M(t) = 



Me 

M4 



t - 1, 



-,1/2 



1/4. 

(3.12) 
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Here we took into account the appropriate initial con- 
dition that M — Mc Sit t — tc- From this formula the 
lifetime of black hole icvap with initial mass Mc can be 
estimated as 



^cvap — 



ti + 



(3.13) 



Note that the lifetime for a PBH hotter than 100 MeV is 
significantly shortened and the above expression cannot 
be applied any longer. 

We derive some typical values of black hole mass and 
temperature. First of all, PBHs just formed at the end of 
the brane high-energy phase have mass Mc I and radius 
rg ~ I. Then I ^ 10^''?4 is needed to lower the temper- 
ature of such PBHs around 100 MeV and hence above 
10^^/4 the preceding estimation of lifetime is valid; it 
gives iovap ~ {l/10'^%)Ho. Thus we find that PBHs evap- 
orating at the present epoch enjoy five-dimensionality 
only with I larger than 10^*^/4. Later on it will also be re- 
vealed that PBHs relevant for our analysis are only those 
with a lifetime ~ to- If ^ falls below 10^*^/4 we cannot see 
the difference from the 4D case. Therefore we only have 
to pay attention to the case I > lO^'^U throughout this 
paper. 



For those black holes with lifetime tr 



^ tn ^ 



^cvap ''0 

13.7 Gyr, their initial mass M* is calculated as 



AC = 3.0 X 10^ 



1031L 



-1/2 



Thus corresponding initial temperature is 



^BH = 57 



I 



1031L 



-1/4 



keV. 



(3.14) 



(3.15) 



With / in the typical range 10 Z4-IO ^4 of our interest, 



57 keV < Tj^H < 32 MeV, 



(3.16) 

hence the consistency of low temperature approxi- 
mation for those PBHs is confirmed. Furthermore, 
since the black hole mass corresponding to Tbh ~ 
100 MeV is M|tbh=ioo MoV 10*(Vl0^^^4)~^ g, the ra- 
tio M|tbh=ioo Mew/M* is much less than unity except I 
is very close to 10^°Z4. Therefore we also confirm that 
the fraction of photon energy emitted after the Hawk- 
ing temperature rises above the cutoff temperature m^r 
is negligible for those long-life PBHs. 

Summarizing the previous sections, our formula to de- 
scribe the time evolution of black hole mass is obtained 
as 

2F/7r 

Mi for t < tc (3.17a) 



Mit) = - 



Me 
M4 



where Mc = M {tc 



t - 

<4 



1/2 



for t > tr 



Ma 



(3.17b) 



D. PBH abundance 

In order to evaluate the fraction of energy density in 
the form of PBHs to that of radiation, we take an usual 
notation 



a* (Mi) 



PPBH^Afi (t) 
Prad(0 



(3.18) 



where ppBH,Mi is the mass density of PBHs with forma- 
tion mass on the order of Mi, and piad is the radiation 
energy density. In this paper we only use ai = at^ , where 
ti is given for Mi by solving Eq. 13.1|l . at is related to ai 
via cosmological scale factor a{t) as 



at{Mi)^ai{Mi 



aft) 



(3.19) 



ai is an important quantity for the reason that it is di- 
rectly connected with the primordial density perturba- 
tion. The upper limit on ai, which we are seeking, is 
written as LIMi. 

Using the above definitions, the number density of 
PBHs with formation mass on the order of Mi is given 
as 



n{Mi) = Qi(Mi)-^ ^ ^(tcvap 



(3.20) 



where tovap is corresponding lifetime and 9(t) is the theta 
function. 

We note that the above number density n{Mi) is a fi- 
nite quantity although it apparently depends on specific 
mass A/i; this is because the mass fraction Q!i(Afi) is a 
quantity for PBHs with mass on the order of Mi. The 
relation between n(Afi) and the differential mass spec- 
trum dnpBn/dAfi can be written as 



n(Afi) 



0.1/3 Mi 



^^dnpBH . . 



where /? is a constant on the order of unity. Therefore 
some constraint on n{Mi) leads to an equivalent con- 
straint on the mass spectrum dnpBn/dMi. 

Throughout this paper, we only consider the number 
density n(Mi) rather than the fundamental mass spec- 
trum. This treatment, although it has been usually em- 
ployed in arguments on PBHs, may cause differences in 
the results due to the uncertainty in the definition of 
n(Mi), however, our conclusions will not be qualitatively 
changed. In the following discussion we often call num- 
ber density n(Mi) "mass spectrum" because those two 
quantities are essentially equivalent. 



E. Diffuse photon 

As a consequence of Hawking radiation, PBHs emit 
particles at every moment and then reduce their masses 
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obeying the formula (|3.17b|) . Particularly we are inter- 
ested in photons emitted between the decoupling time 
tdcc and the present epoch to because they have never 
strongly interacted with matter and hence its original 
spectral shape should be preserved in the observable 
high-energy diffuse photon background. But unfortu- 
nately it is understood that we cannot see its signature 
in the present observations; all we can do is to derive ex- 
pected spectrum and set the upper bound on PBH abun- 
dance. 

As for photons, the relevant differential emission rate 
of Hawking radiation is 



dEdt ~ exp(£;/TBH) - 1 ^ 



(3.22) 



Instead of integrating it, we use an approximated expres- 
sion for the emission rate of energetic photons: 



dN ^ d^N 
E- 



dt dEdt 



(3.23) 



Thus total emission rate of energy density via photons 
with energy E is given as 



dt 



dEdt 



(3.24) 



In the cosmological context, photon energy E emitted at 
time t suffers from redshift. Thus the observed energy of 
photons emitted at time t with energy E is 



a[toj 



(3.25) 



at present. Additionally, the energy density of photons 
is also reduced by a factor of a~*. Therefore resulting 
energy density Uo{Eo) observable now is expressed as fol- 
lows: 

= fdM^'-^r-dtE^"''' 



dMi 



dEdt a4(io) ■ 

(3.26) 



Here the effect of absorption during propagation in the 
extra-galactic field is ignored because it should not cause 
significant deviation in the energy range of our interest. 

Since the energy spectrum of emitted particles at a 
moment shows a sharp peak, it is generally expected that 
the spectrum of Uq has two peaks, or at least edges. One, 
the high-energy peak, is around the temperature of PBHs 
with lifetime to. The other, the low-energy peak, is also 
determined by the temperature of PBHs with lifetime idoc 
but its location is redshifted by a factor of a(idcc)/a(^o)- 
The high-energy peak takes a crucial role in the following 
discussions, but the low-energy peak is less important. 



When we are to compare the calculation with the ob- 
served spectrum, an adequate quantity is the surface 
brightness. 



I{Eo) 



c UojEp) 
An Eo 



(3.27) 



From now 

on, "(diffuse) photon spectrum" is used for this quantity 
throughout this paper. 



IV. ANALYSIS ON SPECTRA 

In this section, we investigate two different kinds of 
spectra, the PBH mass spectrum and the diffuse photon 
spectrum. 

To describe the time evolution of mass spectrum, we 
have to consider two kinds of masses, Mj and Mc- Al- 
though we defined the former as the mass at the ini- 
tial time of individual PBHs, it appears natural to treat 
the latter as an initial value of successive time evolu- 
tion. Therefore in the following sections we basically call 
Mc and Mj "initial mass" and "formation mass," respec- 
tively. 



A. Relation between PBH mass spectrum and 
diffuse photon spectrum 

Since we are mainly interested in PBHs having evapo- 
rated away until now, we consider light PBHs formed in 
the brane high energy phase, where the approximation of 
ti < tc holds weU. From Eqs. (|Oa|) . ljI!T|l . and llOHl . 
the number density of PBHs with formation mass Mi is 
obtained as 



n(Mi)|t 



a; (Mi 
3M| 



Prad(^i) 
Mi 



4tt 



ai(Mi)/i/2r3/2Mr3/2. (4.1) 



We should recall that this quantity n\t^t^ is not a spec- 
trum at a specific time but only a collection of the values 
of number densities evaluated at the formation time of 
each PBH. Thus the above mass spectrum is distorted by 
subsequent accretion; the mass growth ratio for individ- 
ual PBHs differs according to their formation mass Mi. 
Therefore resulting mass Mc after accretion is obtained 
from Eqs. (|Oal) . lEUJ, and (|3.17a|l as 



Mc 



2F/TI 



Mi 



4M; 



Mi) 



F/-n 



Mi. 



Thus in turn the following relation holds: 



Mi = (AMlfl) 



-F/iTT-F) 



^jn/(-n~F)^ 



(4.2) 



(4.3) 
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Combining Eqs. (|4.1|l and H4.3|l . we obtain the initial 
PBH mass spectrum at the end of accretion as^ 



3M 



9/4 



27/2, 



-ai(Afi)/i/«ri5/8Mr9/8 



3M, 



9/4 



(4.4) 



where a new accretion parameter = 9F/ (Stt — 8F) was 
introduced for convenience. Though there remains 
as the argument of a;, it is exphcitly determined by 
using Eq. (|4.3(l . Finally, mass spectrum at any later time 
is given by the following form: 



n(Me) = ^^n(Me)|t=t, 



a{ty 



cq 



X (4M2)^ j1/8+^;-3/8+^^^-9/8-^_ (4 5) 

Particularly we find a relation between the distorted 
number density and that with no accretion as follows: 



n(Afe) = 



Qi(Mi) 
ai{Mc) 



"^^lijPj '^(^^c)|„o accretion. (4.6) 



We realize that the accretion controls how effectively 
other braneworld parameters I and / modify the mass 
spectrum of PBHs. 

Now we evaluate the parameter-dependent factors in 
the above expression. We first recall that variables that 
have uncertainties on the orders of magnitude are ai and 
I. On the other hand / stays at 0(1). Since F runs 
from zero to around unity, J-^ is typically small: < 
T < 9/(87r — 8) « 0.53. Then it appears that for any 
values of constants and / in the round brackets in 
Eq. H4.6|l have little significance on mass spectrum. Thus 
we concentrate on I- and Mc-dependent parts. 

At first we shall fix I and let F run to investigate the 
pure effect of accretion. We see in Eq. (14.51) that JT, ac- 
cretion itself, can modify the spectrum generally softer: 
from cx ai{Mi)M-^-^ to maximally cx ai{Mi)M-'^-'' . This 
reflects the fact that the number of lighter PBHs is 
more enhanced than that of heavier ones by accretion. 
Whether the number density of PBHs with initial mass 



Mc is enhanced or reduced is determined by the propor- 
tionality 



ai(Mi) / I 



(4.7) 



Because typical values of I and Afc which we are now 
interested are in the ranges 1 <C Afc ^ 10^° < I < lO'^^ in 
4D Planck units, the round bracket gives an enhancing 
factor of orders of magnitude. On the other hand we have 
less knowledge of a\. However, we expect a\ does not 
vary so rapidly in Af, hence we never get involved in a\ 
too much and hereafter omit its contribution. Therefore 
it is concluded that in this case the number of PBHs is 
actually enhanced and it directly leads to the gain of a 
resulting diffuse photon spectrum. 

Next we fix F and let I vary. When I varies over orders 
of magnitude, number density also varies by nearly the 
same orders as n cx However, we should remember 
that varying I results in a change of mass range rele- 
vant to the observable diffuse photon spectrum. Thus we 
derive parameter-dependent behavior of a photon spec- 
trum in the following manner. First the absolute values 
of the present photon spectrum depend mainly on the 
total amount of PBHs having already evaporated away 
by now because most of the energy emission from indi- 
vidual PBHs occurs around the end of their lifetimes. 
Here, since the lifetime tevap is proportional to Af^/, 
the initial mass of PBHs with a specific lifetime obeys 
Afc cx Therefore from Eq. H4.5I) the number of 

PBHs which evaporate at some cosmological time obeys 
n cx cki 

(Afi);3/i6+3.;F/2 Y^otcA here that the abso- 

lute value of the resulting photon spectrum is related not 
only to the number density of evaporating PBHs but to 
their mass and temperature. Now we shall concentrate 
on the peak value of the spectrum and investigate its I 
and F dependences. In order to evaluate the peak, we 
here make an approximation that photons relevant to the 
high-energy peak are only those emitted from PBHs just 
evaporating now.'^ Furthermore, we also regard that the 
number of emitted photons per PBH is proportional to 
Afc/TeH- Thus the dependence of the peak value fpcak 
on parameters is obtained as 



^pcak 



ai(Af*)ri/i6+3^/2_ 



(4.8) 



We further derive the relation between fpcak and 
Hawking temperature Tg^, which describes a trajectory 
of peaks when I is varied but F is fixed. This representa- 
tion will turn out to be greatly useful in later analyses. 



^ Because n is an integrated quantity, Jacobian only contributes 
as a factor of order unity in the final form and we omit it for 
simplicity. From Eq. i3.21i . one can see that this factor takes 
the form of tt / {n — F) . 



^ The intensity o f th e emitted particle shows a peak at ~ 5 X Tbh 
in the 4D case [Tot . 



9 



Substituting I (x Tgjj , the resulting relation is given as 
a,{M*)El, (4.9) 



^pcak 



where we changed the variable from Tgjj to E'o, and de- 
fined p = 1/4 — QT. The region for the parameter Eq 
is limited to 10 keV < Eq < 100 MeV as in Eq. l|TTB)l . 
We find that there is a specific value F « 0.11 that sets 
the above index to be zero. However, in our framework, 
it will appear that a more significant value is F k, 0.52 
giving the index p k, —1.1, which is actually the index 
of the power-law fitting of the diffuse 7-ray background 
observed by EGRET 



10'" 
10" 
•T 10*^ 

l 10^ 



10" 
10-^ 
10-" 
10"* 
10"* 



10" 



-10 



10" 



'T 1 •— TT 1 •— TT 1 .— ry- 



Observed diffuse photon background «£o 



f>0.5 




10^ 10" 10" 
£0 [keV] 



B. Comparing with observations 

From now on, we compare the predicted spectra with 
the observed diffuse photon background and then obtain 
the upper limits on the mass fraction aj by requiring that 
those spectra do not contradict the observation. LIMiS 
are determined for each parameter set. 

First of all, the observed high-energy diffuse photon 
background spectrum in the region 100 keV - 1 GeV 
obeys approximately the power-law with index —1 to —2 
I33L 0, [3^ l36|. On the other hand it was shown by 
Clancy et al. p33 | that a photon spectrum emitted by 
PBHs asymptotically approaches ^ E^^^ in the high en- 
ergy tail, and has typically more hard shape than Eq^ 
in the energy range lower than the peak as long as ai 
does not vary so quickly with Mi. Therefore under the 
same assumption it is understood that the spectral peak 
generally gives an upper bound on ai. Keeping this fact 
in mind, we further simplify so that ai apparently does 
not depend on Mi in the following analysis. 

In this case the peak trajectory Eq. (|4.9|l is fixed except 
with respect to the "normalization factor" ai. Here, con- 
sider the case in which the trajectory and the observed 
background spectrum have an intersection at some large 
tti. Apparently the / corresponding with the intersection 
bounds the size of the extra dimension. At this point we 
realize that it is possible to classify the situations into 
two distinctive classes; one is such that the gradient of 
the peak trajectory p < —1, i.e., F > 0.5, and the other 
is p> —1, F < 0.5. The scheme is shown in Fig. 

In the former case ("F > 0.5"), the peak trajectory 
is "softer" than that of the observed background spec- 
trum. We should recall that higher energy corresponds 
with smaller I. Thus once the value of / at the inter- 
section is given, it sets an upper bound on /; the size of 
the extra dimension I is allowed only if it is smaller than 
this bound so that the peak is hidden in the photon back- 
ground. Saying with the aid of Fig.^ the peak locates in 
the side of (a) rather than (b). On the other hand, in the 
latter case ("F < 0.5"), the trajectory is "harder" than 
the background. It can be even an increasing function 
of Eq. Somewhat surprisingly, this time the intersection 



FIG. 1: A schematic view of the relations between the diffuse 
photon background and two kinds of trajectories. The peak 
energy is limited in the region shown in Eq. 13.161 . thus the 
peak can move between either pair of ends. In each case Oi's 
are set to be large enough to make trajectories and the photon 
background cross. 

gives a lower bound on /. Namely, a rather large extra 
dimension is required, i.e., (c) is allowed but (d) is re- 
jected. However we again should note that the above 
discussions are valid only when the trajectory intersects 
the background. Since Eq moves in a limited range, a 
much smaller ai prevents them from crossing. 

For further quantitative investigation, we consider 
the transition between five dimensions and four dimen- 
sions. To begin with, we determine when and how four- 
dimensionality is recovered in our framework in the fol- 
lowing manner. Although we naturally expect that ev- 
erything should recover four-dimensionality in the small I 
limit, when the recovery does occur is not trivial because 
it depends on what quantity we treat. In our framework, 
the observable quantity is the shape of the photon spec- 
trum emitted from PBHs in the past. There are espe- 
cially two elements, the absolute value and the location 
of the peak. The former depends on both the size of the 
extra dimension / and the accretion efhciency F, while 
the latter is determined by / alone. The peak energy 
leaves the 4D value only when the Hawking temperature 
of PBHs evaporating now differs from the 4D one. The 
difference in Hawking temperature comes from the in- 
equality of the Schwarzschild radius rs of the PBH vs 
the AdS curvature radius If a PBH has sufficiently 
smaller rs than /, it is well described as a 5D object and 
its Hawking temperature deviates from the 4D case ac- 
cording to Eqs. H2.15|l or H3.15|l . Therefore the change 
of peak location can occur when PBHs evaporating now 
have radii smaller than /. However, in the opposite case 
where they have a larger radius than /, their tempera- 
ture and the evolution are the same as in the 4D case. 
Although a spectral shape other than the peak can be 
modified a little by the 5D effect, properties of the peak, 
such as the location and its absolute value, are indistin- 
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guishable from 4D ones. Thus we can set the boundary 
between 5D and 4D on the case when the size of the 
Schwarzschild radius of newly evaporating PBHs agree 
with the size of the extra dimension: rg = I, which cor- 
responds to I « 5.1 X 10^^/4. With an additional fact 
that PBHs formed at t ^ ic have rs ^ I, this condi- 
tion is equivalent to the one that the 5D accretion effect 
becomes significant, as shown in Sec. IIIII 

By the way we have already found that there is some 
region in parameter space with low accretion efficiency 
(F < 0.5) in which a smaller extra dimension is disfa- 
vored. Of course the existence of such a strange region is 
entirely due to the existence of the extra dimension, but 
from the above arguments, it was revealed that there is 
some critical size of I in our framework below which the 
effect of the extra dimension is hidden. Therefore we con- 
clude that the favorable I with low accretion efficiency is 
at least larger than the critical size 5.1 x 10^^^4. 

We can confirm the above implications by taking up- 
per bounds on ai into account specifically. Let us first 
consider the case of F < 0.5, in which the peak trajec- 
tory obeys /peak Eq with the index p > —1. The 
severest bound on ai is set when the peak locates at the 
high energy end of its trajectory, i.e., ~ 100 MeV, and 
comes in touch with the observed background spectrum 
there. Then the effective 4D condition should hold and 
the corresponding upper bound on a-^ is exactly the same 
as that of the 4D case, 10"^'^ On the other hand, if 
the Universe has an extra dimension with a size exceeding 
the critical value, the peak moves to a lower energy re- 
gion sufficiently hidden in the background spectrum and 
hence an ai larger than 10^^^ becomes allowed. As a 
consequence ai larger than lO"^'' requires / to be larger 
than 5.1 x 10^^/4. Note that the above conclusion makes 
sense only when the conditions ai > 10~^^ and F < 0.5 
are simultaneously realized. The former is forbidden in 
the pure 4D universe but the latter means that the PBH 
mass spectrum is almost the same as that of the 4D case, 
where no accretion works. Therefore the existence of a 
large extra dimension, which can decelerate the evolution 
of PBH, is needed to avoid conflict. On the contrary no 
lower bound on I is set when ai is smaller than lO"^''. 

Finally, we mention the case F > 0.5. Like the 
F < 0.5 case, the peak trajectory intersects the observed 
background spectrum at the highest energy end when 
ai = 10~^^. However, this time the trajectory entirely 
locates above the background spectrum except at the in- 
tersection. Thus in order to make the extra dimension 
larger than 5.1 x 10^^/4 and let the peak move to a lower 
energy region ensuring observed results, the trajectory 
must fall down by setting ai lower than 10~^^. Conse- 
quently ai > 10~^^ is never allowed and I larger than 
5.1 X 10^^/4 gives a rather stringent upper bound on ai 
than in the 4D case. Of course there is no deviation from 
the 4D case if / is below the critical size. 



V. NUMERICAL RESULTS 

Here we confirm and refine the arguments in the previ- 
ous section by virtue of numerical calculation. Relevant 
formulas are Eqs. and (|T7r|) with Eq. At 

first, we study how the spectral shape of the PBH pho- 
ton depends on each kind of parameter, and then obser- 
vationally allowed regions for braneworld parameters and 
their implications are shown. 

In the following calculations, we assume that ai is con- 
stant in All for simplicity. We finish calculating the evap- 
oration process of the individual PBH when its Hawking 
temperature reaches ~ 135 MeV. 

Observational data being used in this section are listed 
in Table m We should note that there is relatively large 
error in the range 400 keV-30 MeV among these obser- 
vations. There remains uncertainty in our results coming 
from it. 



TABLE I: Observational data used for patches of energy 
regions and their references. 



3 keV 
80 keV 
450 keV 
- 1 MeV 



60 keV 
450 keV 
~ 1 MeV 



HEAO-1 A2 [33j 
HEAO-1 A4 [34] 
SMM'-tSSj 
EGRET [36] 



"There seems to be no appropriate data in 60-80 keV, but this 
energy range is not a concern in our discussions. We just draw the 
interpolated spectrum in this region. 

''Since there are relatively large uncertainties in this energy region, 
we use only the index of SMM and connect the values of HEAO- 
1 A4 and EGRET. This method may cause an error of an order 
of magnitude in our conclu sion . Note that such a spectrum is 
consistent with COMPTEL Isl. 



A. Accretion efficiency F 

We first discuss the effect of accretion on PBHs in the 
brane high energy phase. Changes in the spectra caused 
by varying F are shown in Fig. Accretion efficiency 
F is set to 0, 0.5, and 1.0, while other parameters are 
fixed to ai = 10~^^, I — lO^^l^, and / = 1. Because 
other parameters I and / are fixed, the initial PBH mass 
spectrum and subsequent evolution of individual PBHs 
are both common. In addition, the mass range relevant 
to the resulting diffuse photon spectrum, namely of PBHs 
with lifetime fdcc £ ^cvap < ^Oi is never changed. 

The spectra show two peaks as mentioned in Sec. IIIII 
The high-energy peak is around (a few) x 100 keV, which 



* Note that hereafter we set the value of the initial mass fraction cti 
to be lO"'^^. Those spectra crossing the observed diffuse photon 
background are observationally rejected for this standard oj, e.g., 
the F = 1.0 case in Fig.|21 



11 



10'" 
10* 
10* 
10^ 
10^ 
10° 
10-^ 
10-" 
10"' 
10"* 
10-'° 





11,1 I 1 1 I ' ' 1 ' 

ai=10--', /=10^'Z4,/=1 - 


f=1.0 








. . 1 . . . 1 . . . 1 . V^, 1 , , "."i-. , , , 



10" 



10' 



10' 



10^ 
Eo [keV] 



10" 



10" 



10° 



FIG. 2: Spectra of surface brightness / for F — (long- 
dashed), 0.5 (short-dashed), and 1.0 (dotted), respectively. 
Other parameters are fixed: ai — lO"^'^, I = 10"^^Z4, and 
/ = 1. The solid line indicates an observed diffuse photon 
background. Spectra with accretion {F = 1.0 and 0.5) are 
almost uniformly gained according to Eq. (14.711 relative to that 
with no accretion {F = 0). A slight difference of gradients 
in the low energy region between these spectra is seen; it 
comes from the fact that the number of lighter mass PBHs 
is relatively more enhanced by accretion than that of heavier 
PBHs. 



corresponds to the Hawking temperature of a PBH with 
hfetime to. See Eq. (|TT^ and substitute / = 10^^;4. The 
enhancement of absolute value in this case comes from 
the number growth of PBHs as in Eq. (|4.7|l . Specifically, 
making use of Eq. H3.14|l with I — 10'^^^4, we obtain fac- 



tors of 



I 



= m 



8.9 



10' 



3.6 



for F = 1.0 



for F = 0.5, 



(5.1a) 
(5.1b) 



which explain calculated spectra well. 

The change of gradient in the lower energy region is 
due to the distortion of the initial mass spectrum. In 
short because of relatively large enhancement of smaller 
PBHs. 



B. Bulk AdS curvature radius I 

Next we consider the dependence on bulk AdS curva- 
ture radius I, one of the most fundamental quantities in 
the braneworld scenario. 

At first glance, one finds two basic facts, that (i) the 
Hawking temperature of a black hole on the brane is low- 
ered as oc l^^/^ in contrast to the 4D case, and that (ii) 
the brane high energy phase (t < 1/2) becomes longer 
for larger I. Thus in general the behavior is rather com- 
plicated. However, we can neglect the second fact when 
accretion does not work. Hence we first concentrate on 



cases in which accretion never works {F = 0), then we 
consider the composite cases {F = 1.0 and 0.5) later. 
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FIG. 3: Spectra of surface brightness I for I = 10^^/4 (long- 
dashed), 10^^/4 (short-dashed), and 10'^^Z4 (dotted), respec- 
tively. The other parameters are fixed: Oi = 10~^^, F = 0, 
and / = 1. A 4D spectrum (thinsolid line) is also shown. The 
high-energy spectral peaks move from high energy (right) to 
low energy (left) as the extra dimension becomes large. Ac- 
tually, the peaks draw the trajectory /peak tx Eq'^^ according 
to Eq. 14.91 . As seen in this figure, it is expected that the 
upper bound on Oi for this parameter set of F and / is given 
in the high energy region, and thus in turn large I is favored. 
Convergence on the 4D spectrum as I approaches 10^" U is 
also shown. 



Figure 131 shows dependence of spectral shape on AdS 
curvature radius I in the case of no accretion. I is taken to 
be 10^1/4, 10^6/4, and W^Hr, f and F are fixed to unity 
and zero, respectively. The 4D spectrum is also shown in 
the figure. 

At first we restrict ourselves to the three 5D spectra. 
Since there is no accretion, mass spectra of PBHs in each 
case are identical. They are expressed as Eq. I)4.5|l with 
T = 0. The only difference between the three spectra 
comes from the Hawking temperature, equivalently the 
lifetime of PBH. For larger I, the Hawking temperature of 
PBHs is lower. Thus the number of PBHs contributing to 
the peak is enhanced, but on the other hand the number 
of photons emitted from individual PBHs is decreased. 
The two opposite effects almost compensate in this case, 
see Eq. (|4.8|l and set = 0. As a consequence, the 
peaks draw an almost horizontal trajectory /peak oc Eq''^^ 
according to Eq. (|4.9|l . 

Next we pay attention to the 4D spectrum. We can 
confirm the consistency of our calculation with the 4D 
case in Fig. |31 from the converging behavior of spectra 
upon the 4D spectrum as I becomes small. As mentioned 
in the previous section, the transition from 5D to 4D 
occurs when the 5D approximation for the PBHs with 
lifetime to breaks down. The critical size of / is w 5.1 x 
10^^^4. The properties of PBHs with lifetime tovap = 
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under such a condition, 

Tj^H = 38 MeV and M* = 1.3 x 10^^ g, (5.2) 

are both consistent with known 4D values. In addi- 
tion, under this 4D hmit, the brane high energy phase 
is shortened and number growth of PBHs hardly occurs. 
As a consequence 5D spectra succeedingly reproduce the 
shape of 4D spectrum with respect to both the location 
of its peak and absolute value under the 4D limit. 

Note that the 4D spectrum sets the upper limit on 
the mass fraction ai around ^ 10~^^ We already 

argued in the previous section that the value 10~^^ takes 
a crucial role even in the braneworld. We later confirm 
it numerically. 

By the way, now we are ready to understand the gen- 
eral behavior of spectral shape for varying I together with 
the existence of accretion. See Fig. 0] 

The locations of peak energy are explained in the same 
manner as in the case of no accretion. In addition, the 
dependence of absolute values on parameter sets is ex- 
plained by Eq. H4.8|l as 

/peak = Z"-^' for F = 1.0 (5.3a) 
= fori^ = 0.5, (5.3b) 

and the trajectories of their peaks are 

/peak - ' for F = 1.0 (5.4a) 
= E^''-^ for F = 0.5. (5.4b) 

We can see that in the case F = 1.0 an upper bound 
on I is given at the intersection of the above trajectory 
and the observed background spectrum. To the contrary, 
in the case = 0, a lower bound is obtained from the 
intersection. On the other hand the spectral index for 
F = 0.5 is nearly identical to that of the observed back- 
ground. Thus we confirm the mention in the previous 
section that the bound on I switches its character as an 
upper bound to that as a lower bound when accretion 
efficiency falls down across 50%. 

C. Uncertainty of formation process 

The horizon mass fraction parameter / takes a less 
important role in the results. It only contributes as an 
overall factor for the number density of PBHs in the form 
n oc as shown in Eq. H4.5|l . From a physical point 

of view, / has a parameter range at most within a few 
orders of magnitude. Therefore we skip showing the spec- 
tral shapes for varying /; they are not so interesting. 

D. Allowed regions: The final result 

As shown in this section, we have obtained photon 
spectra for sufficiently many parameter sets. Here we 



apply those sample spectra to determine how large the 
initial mass fraction ai can be. 

To achieve it, we adjust the calculated spectra to come 
in contact with the observational data, then we can de- 
termine the value of the upper limit on ai, LIMj, by mea- 
suring overall factors of adjusted spectra. The method is 
shown in Fig. \^ 

It is confirmed in Fig. \S[ that the most crucial part of 
each spectrum giving the upper limit on a; is nothing but 
the high-energy peak. 

Collecting values of LIMi for _F's, we can show the in- 
dependence of LIMj as in Fig. |H1 where we selected two 
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FIG. 4: Spectra of surface brightness / for I = 10^^/4 (long- 
dashed), 10^^/4 (short-dashed), and 10'^^/4 (dotted), respec- 
tively, where the other parameters are fixed: ai — 10~^^, 
/ = 1, (a) _F = 1.0, and (b) 0.5. In both figures, the spectral 
peaks or edges move from the high energy region (right) to the 
low energy region (left) as the extra dimension becomes large. 
The peaks in two figures draw trajectories /peak oc E^^'^ [(a) 
F = 1.0] and Eg'^-° [(b) F = 0.5] according to Eq. (ll^J. By 
considering the gradients of trajectories, we find that upper 
bounds on I are given for _F = 1.0 from the value of / corre- 
sponding with the possible intersection of the observation and 
trajectory; peaks must locate in the right of the intersection. 
In particular, ai < 10~^^ is required to realize i > 5.1 x 10^^/4. 
Note that in the F = 0.5 case, the gradient of the peak trajec- 
tory and that of the observed spectrum are almost the same. 
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FIG. 5: Illustration of how to obtain LIMi's. The six spectra 
in this figure are those that appeared in Fig. |3 and Fig.|lja), 
the F = 1.0 case, but this time their absolute values are 
adjusted to be in contact with the observed spectrum, i.e., 
Oi 7^ 10~^^. These qi's give upper bounds on the fraction of 
PBHs for each parameter set. 




FIG. 6: Illustration of how LIMi's depend on _F's in the two 
cases / = 1 (solid) and / = 0.001 (dashed). For each /, three 
sizes of extra dimension are taken: / = W^^U, W^^k, and 
10^^ Z4 from above to below in the figure. The differences of 
qualitative behavior between the lines are coming from the 
sensitiveness to accretion efficiency. For a larger extra dimen- 
sion, the duration of the brane high energy phase is length- 
ened and hence the effect of accretion is stressed. To the 
contrary, as the size of the extra dimension becomes small ac- 
cretion ceases to work for any values of F, and LIMi becomes 
less sensitive to _F. LIMi is almost independent of F under 
the limit / 10^° h- 



/'s, 1 and 0.001. In this figure we find that LIMj's fall 
as F's become large for almost all the parameter sets 
displayed here, while one set I = 10^^/4 and / = 0.001 
shows subtle behavior. It should come from a detailed 
spectral shape around the peak, thus beyond our approx- 
imated analysis. When accretion in the brane high en- 



ergy phase works efficiently, the number of PBHs grows 
along with increasing F, and it leads directly to a gained 
peak. Therefore the upper limit on aj is lowered in order 
to match observations. On the other hand as I becomes 
small, accretion docs not have an influence on the result- 
ing spectra and then LIMj becomes independent of F. 
Again we confirm here the importance of 50% accretion 
efficiency in Fig. El which can be read from the behav- 
iors of lines; the lines all pass through around the point 
F = 0.5 and LIM; = lO'^^. 

The constraint LIMi = lO"^'^ obtained by Clancy et 
al. [sj corresponds with the point {F = 0, LIMi ~ 3.9 x 
10~^^) on a line with I = lO^^U and / = 1. Apparently 
the constraint is relaxed by some numerical factor for this 
case. 

At last, we obtain the final result. Until now we have 
taken a stance on searching for upper limits on the mass 
fraction of PBHs. However, once the actual value of ai is 
obtained or constrained in some way, one can determine 
whether some parameter combination {l,F,f) is allowed 
or not by comparing it with LIMi of the set. Namely the 
combination will be observationally rejected if cti > LIMi. 
Thus we obtain a tool for limiting the braneworld param- 
eters through the future knowledge on ai. We indicate 
the diagram with boundaries for ten ai cases in Fig. [7| 
Each boundary separates the parameter space into two 
regions; the left side is allowed and the opposite is ex- 
cluded. 
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FIG. 7: Showing of boundaries for Oi = 10"^^-10"^^ from 
left to right. The remaining parameter / is here fixed to unity. 
The critical size of extra dimension at which transition from 
5D to 4D occurs is also shown (dotted line). The experimen- 
tally excluded region I > 10^^ k is shaded. The boundaries 
split the parameter plane into two regions, respectively, the 
left side is allowed and the right side is excluded. Once PBHs' 
mass fraction Oi is determined in some way, this diagram in- 
dicates what region in the space of braneworld parameters is 
allowed. Further information obtainable from this diagram is 
argued in the text. 

Let us investigate this diagram further. To begin with, 
we should confirm the consistency of this diagram with 
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the 4D case, in which LIMj ~ 10~^ was obtained by 
Green and Liddle H^. As was confirmed in Fig. El we 
expect that our results under the 4D limit should be inde- 
pendent of F. Actually we see such a signature in Fig.[7| 
from the fact that all boundaries go to either the F = 
axis or = 1 before they reach Z = 5.1 x 10^^/4. The 
above qualitative behaviors denote that under the 4D 
limit I < 5.1 X 10^^/4 all the sets of braneworld parame- 
ters with ai > 10~^^ are rejected, while with ai < 10^^'' 
are all allowed. Therefore we find that our result under 
the 4D limit is well consistent with the pure 4D result 
LIMi ~ 10-27. 

Furthermore, there is a critical boundary between ai = 
10^2^ and 10^27, which separates the whole region into 
two distinctive patches. Note that we can identify this 
critical boundary with a vertical line F « 0.5. In the 
left side of the critical boundary, that is ai > 10^^^, all 
boundaries hit F — 0, while for Ui < 10"^'', all hit F = 1. 
Therefore we find that only the region Z > 5.1 x 10^^/4, 
F < 0.5 is allowed for cti's higher than the critical value. 
It means that if large a^, which is not allowed in the 
4D case, is actually realized, large extra dimension is 
required and I has a lower bound > 5.1 x lO^^l^. On the 
other hand the boundary with ai = 10-^^ asymptotically 
approaches to I — 5.1 x 10^^Z4. This means that an extra 
dimension with / larger than 5.1 x 10"'^^^4 requires ai to be 
smaller than 10^^^. Hence a large extra dimension gives 
more stringent upper bound on ai than the 4D case. Of 
course modestly large ai in turn gives a rather severe 
bound on the size of extra dimension. For example, if 
ai = 10^27 with F = 1.0 is realized, I is bounded below 
lO^^U, which is about 11 orders of magnitude severer 
than the present bound obtained from measurement of 
the short-range gravitational force. However, one should 
always keep in mind that if there is merely a smaller 
initial PBH fraction than 10^^^, our framework based on 
the diffuse high-energy photon background has limited 
significance on the braneworld parameters. Even in the 
most efficient accretion case F = 1.0, the small fraction 
10-31 only gives I < 10'^°Z4, which is almost identical 
with the known upper bound. Moreover, F < 0.5 gives 
nothing new if ai < lO^^^". 

VI. CONCLUSIONS AND DISCUSSIONS 

In this paper, at first we started with the mass spec- 
trum of PBHs at the moments of their formation. Then 
we derived the resulting distorted spectrum after accre- 
tion in the brane high energy phase. We found that the 
PBH number density with specific initial mass Mc is en- 
hanced relative to its no-accretion value as 

n=^7^f4M|^') n|„o accretion, (6.1) 

ai(Mc) V McJ 

where 

Mi = {AMjfiy^'^''-^^ M;/(--^). (6.2) 



Using the result on mass spectrum distortion, we also 
derived the behavior of the peak value of the diffuse pho- 
ton spectrum. Under the assumption that ai(Mi) is con- 
stant, the peak value depends on Z, or equivalently on 
Tgjj as 

/pcakOc/-l/16+3^/Vr*H'/*-*^^. (6.3) 

The latter proportionality is rather useful when we com- 
pare the peak trajectory with the observed power-law dif- 
fuse photon background with index « — 1. The exponent 
of the peak trajectory p = 1/4 — gives a value smaller 
than —1 when F > 50%, while it gives a value larger than 
— 1 when F < 50%. Thus there is a critical accretion ef- 
ficiency around 50%, a fairly neutral value. When F 
exceeds 50%, smaller extra dimension and smaller mass 
fraction ai, typically below the 4D upper bound 10-^'^, 
are favored; it is a reasonable situation for the braneworld 
scenario. To the contrary, if F falls below 50%, large 
extra dimension is favored. Furthermore, if actually the 
mass fraction has a value larger than the 4D upper bound 
10-27, the lower bound on the size of extra dimension I 
is even set around 5.1 x 10^^/4. This critical value comes 
from the condition that the Schwarzschild radius of PBHs 
evaporating now is equal to the size of extra dimension. 
We should note that the braneworld signature vanishes 
regardless of the accretion efficiency provided the extra 
dimension is smaller than the above critical size. Nat- 
urally in the case ai ^ 10^^^ ^}jg existence of an extra 
dimension again has nothing to do with our framework. 

In this paper we also demonstrated numerical calcu- 
lations for further preciseness of the above arguments. 
In the calculation we achieved mainly two improvements 
as follows. First, we used the newest observational cos- 
mological parameters obtained by WMAP, which claims 
flA.o = 0.73, n^fih"^ = 0.135, h = 0.71, and to = 
13.7 Gyr. Because the observable diffuse photons are 
emitted by PBHs which evaporated after the decoupling 
time, A-domination in the present epoch could affect the 
resulting spectrum via the accelerated expansion devi- 
ated from the Einstein-de Sitter universe (EdS). How- 
ever, we found that once a correct assumption on the 
present age of the Universe to « 13.7 Gyr is made, the 
difference of scale factor after the decoupling time is at 
most only (a few) x 10% between the two cosmological 
models. As a consequence, the deviation of resulting pho- 
ton spectra was within an order of magnitude. Second, 
we precisely took into account the photon emission pro- 
cess and mass decrease of black holes at every moment. 
It also could modify the result because the redshift is de- 
termined by how long the emitted photons have survived. 
In spite of this effect, the resulting spectrum again does 
not greatly differ from calculations under an approxima- 
tion of instant evaporation. 

After all we obtained numerically precise spectra with 
an assumption of constant ai in sufficiently wide regions 
for each braneworld parameter and compared them with 
relevant observational data. Some of the calculated spec- 
tra were shown in Figs. The behavior of the peak 



15 



value and the constraints coming from them were both 
well explained by the preceding analysis. We then ob- 
tained upper bounds on mass fraction, LIMi's for each 
parameter set {I, F, /) by the method shown in Fig. |S1 
LIMj was sketched in Fig. |S1 Our most impressive re- 
sult is the allowed region in parameter space as shown in 
Fig. The characters of this diagram were also consis- 
tent with our analysis. 

Recently some authors argued a possible modifi- 
cation of the Hawking radiation in the context of 
AdS/conformal field theory (CFT) correspondence in the 
RS2 braneworld |l3)ll3ll3- They predicted enormously 
rapid decay of astrophysical black holes due to emission 
of CFT modes, which is interpreted as the Kaluza-Klein 
(KK) graviton in the bulk. However, the applicable range 
of this effect is rather limited to those PBHs larger than 
AdS curvature radius I. Hence, although the transition 
process of a black hole from 4D to 5D is still not clear, 
large PBHs are to shed their masses and accumulate in 
the mass spectrum around rg ^ Z. They may enhance 
the number of PBHs which can contribute to a diffuse 
photon spectrum and result in a more restrictive upper 
bound on PBH abundance. However, by the following 
estimation it appears that little change is caused in our 
results. We here assume that all PBHs larger than I lose 
their masses shortly after t = and eventually become 
5D black holes with ^ M I, whereas the evapora- 
tion process of those PBHs smaller than I is not affected 
at all. This results in the number growth of PBHs of 
size I, but since in general the realistic number density 
n{M) is expected to be steeper than M"^ though there 
is unknown mass dependence in ai, the enhancement of 
n(A/)|jv,/~i is at most of some numerical factor on the 
order of unity. Moreover, since PBHs with rg ^ I have 
lifetime ievap ~ (VlO^°^4)'^^o, they have little room to 
contribute to diffuse photon spectrum with / in the range 
we concerned {10'^°U < I < lO^^U). 

Although out of our interest in this paper, we point out 
a possibly meaningful case in the framework of AdS/CFT 
mHlia where I < lO'^^U. Under the same assump- 



tions above, in this case we cannot only observe just evap- 
orating 5D PBHs but also so doing standard 4D ones; 
they must have entirely decayed before the present epoch. 
Therefore it becomes more difficult to detect diffuse pho- 
tons emitted by PBHs. Another special case is with crit- 
ical size of extra dimension, / ~ 10^*^/4. This time the 
effect of accumulation around M ^ I can take a signif- 
icant role because they should be evaporating just now. 
However, at the same time they have a rather high initial 
temperature ~ (a few) x 10 MeV and then a fraction of 
their mass energy poured into photon energy is limited. 
It will be necessary to estimate both the compensating 
two effects. 

At last, we mention the simplifications in this paper. 
First, we did not take significantly braneworld inflation 
models into consideration. Though we sometimes took an 
assumption that ai is constant in Mi, the mass spectrum 
may have a rather complicated shape in reality. However, 
if modification is not so radical, it should not affect our 
argument on braneworld constraints because the most 
crucial role is taken by those PBHs which evaporate just 
now. Second, we ignored the latest time behavior of evap- 
orating black holes. It may cause an interesting change in 
spectra, but almost all energy has been already released 
before that. Hence it is not so important for the whole 
discussion above, either. 
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